Entropic elasticity of DNA with a permanent kink 
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Abstract 

Many proteins interact with and deform double-stranded DNA in cells. Single-molecule 
experiments have studied the elasticity of DNA with helix-deforming proteins, including proteins 
that bend DNA. These experiments increase the need for theories of DNA elasticity which 
include helix-deforming proteins. Previous theoretical work on bent DNA has examined a long 
DNA molecule with many nonspecifically binding proteins. However, recent experiments used 
relatively short DNA molecules with a single, well-defined bend site. Here we develop a simple, 
theoretical description of the effect of a single bend. We then include the description of the bend 
in the finite worm like chain model (FWLC) of short DNA molecules attached to beads. We 
predict how the DNA force-extension relation changes due to formation of a single permanent 
kink, at all values of the applied stretching force. Our predictions show that high-resolution 
single-molecule experiments could determine the bend angle induced upon protein binding. 

Keywords: DNA elasticity, force-extension measurements, helix-deforming proteins, transcrip- 
tion factors, bent DNA, theory. 

1 Introduction 

In cells, many different interactions between DNA and proteins occur, processes which are essential 
to gene expression, genome replication, and cellular DNA management. One major class of proteins 
interacts with DNA and mechanically deforms the double helix by wrapping, looping, twisting, or 
bending DNA (0, 0) • Examples include DNA-packaging proteins and transcription factors which 
regulate gene expression. The mechanical deformation of the DNA may be important for gene 
expression: it has been suggested that DNA deformation by transcription factors may help other 
proteins bind to the DNA and initiate transcription. 

The deformation of DNA by proteins can be detected in single-molecule force microscopy. In 
this experimental method, force is applied to individual DNA molecules and the DNA end-to-end 
extension is measured (figure . Single-molecule force microscopy has been used to detect the 
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deformation of DNA caused by protein binding (0, 0, 0, 0, El @). In these experiments the DNA 
end-to-end extension changes when a deformation- inducing protein binds. Varying the applied 
force allows one to probe the deformation and better understand the details of the protein-DNA 
interaction. 

In this paper we focus on proteins that bend the DNA backbone and develop theoretical pre- 
dictions of the force-extension behavior of bent DNA. Our description is based on the worm-like 
chain theory (WLC) @,0,E3)- 

The WLC predicts the average end-to-end extension z of a semi- 
flexible polymer, given the force F applied to the ends of the chain and the values of two constant 
parameters (the contour length L and the persistence length A). However, DNA elastic behavior 
is altered by backbone-deforming proteins, an effect that is not included in the traditional WLC. 
Extended theories have been developed which combine the WLC treatment of DNA elasticity with 
local bends. Rivetti et al. addressed the case of zero applied force |l2h. while Yan and Marko 
have described the changes in the force-extension behavior of a long polymer to which many kink- 
inducing proteins can bind nonspecifically (jlil). Similarly, Popov and Tkachenko studied the effects 
of a large number of reversible kinks |l4h. Metzler et al. studied loop s formed by slip-rings (fl^l ) , 
and Kulic et al. studied the high- force limit of a kinked polymer (|16l ). 

Previous theoretical work has focused on large numbers of reversible kinks or the limit of low 
or high applied force. However, recent single-molecule experiments have examined relatively short 
DNA molecules with a single specific kink site, over a range of applied force (0). Therefore a theory 
is needed which applies to (i) one kink site and (ii) a polymer of finite contour length {L/A ~ 1 — 10). 
Recently, we introduced a modified solution of the WLC applicable to polymers of this length, and 
demonstrated that a pply ing; the traditional WLC solution to molecules with L/A ~ 1 — 10 can lead 
to significant errors (|17h. Our finite worm-like chain solution (FWLC) includes both finite-length 
effects, often neglected in WLC calculations, and the effect of the rotational fluctuations of a bead 
attached to the end of the chain. 



[Figure 1 about here.] 



Here we formulate a theoretical description of a single kink induced by a protein, and extend the 
FWLC treatment to include such local distortions. Our theory has a simple analytical formulation 
for the case of a force-independent bend angle, i.e., a rigid protein-DNA complex. Our predictions 
are relevant to experiments like those of Dixit et al. (0), which detect with high resolution a single 
bend induced in a relatively short DNA molecule. Although we will primarily focus on the case 
of a single bend angle, our method can also describe a kink which takes on different angles with 
different probabilities. This model could be relevant to a binding protein that can fluctuate between 
different binding conformations with different kink angles (flih. 



2 Theory 

[Figure 2 about here. 



2.1 FWLC theory of unkinked DNA 

The classic WLC model (10) and the FWLC theory (17), which includes finite-length effects, de- 
scribe an inextensible polymer with isotropic bending rigidity. The bending rigidity is characterized 
by the persistence length, A, the length scale over which thermal fluctuations randomize the chain 
orientation. We assume that the twist is unconstrained and can be neglected (as is the case, for 
example, in optical tweezer experiments). 
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The chain energy function includes terms which represent the bending energy and the work 
done by the applied force: 
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E = J ' fa U~ /i " t J' 

where E is the energy divided by the thermal energy ksT, t = L/A, s denotes arc length divided by 
the persistence length A, and all other lengths are similarly measured in units of A. The quantity 
/ is the force multiplied by A/kgT, and we assume the force is applied in the z direction. The 
total extension of the chain is z = J ds z • t. The curvature k can be defined in terms of arc-length 
derivatives of the chain coordinate (figure dJl). If the chain conformation is described by a space 

curve r(s) and the unit vector tangent to the chain is t(s), then k 
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If the ends of the chain are held at fixed orientations, we have 

Z(i f ,£;ii,0) = J Dt exp - J ds {^{d s tf - /z • tj , (2) 

where the integral in Dt is over all possible paths between the two endpoints of the chain with the 
specified orientations. The partition function can be interpreted as a propagator which connects the 
probability distribution for the tangent vector at point s, ^(t, s) to the same probability distribution 
at point s': 

ifj(t,s) = J dt' Z(t,s;i',s') ip(i',s'). (3) 

From this relation, one can derive a Schrodinger-like equation, which describes the s evolution of 
V> (0): 

+ (4, 

Here V 2 is the two-dimensional Laplacian on the surface of the unit sphere and cos# = z • t. 

For relatively short DNA molecules (£ ~ 1 — 10), the boundary conditions at the ends of the chain 
(0,0) and bead rotational fluctuations become important. The boundary conditions are specified 
by two probability density functions, ip(t, s = 0) and ^(t, s = L). The boundary conditions modify 
the force-extension relation, and enter the full partition function matrix element via 

Z tot = Jdii dt f 4>(ti, 0)Z(t;,0;t/, L)V(t/,L). (5) 

Rotational fluctuations of the bead(s) attached to the end of the DNA complicate the analysis of 
experiments. What is observed and controlled is not the endpoint of the (invisible) DNA chain, 
but rather the bead's center. The relation between these distinct points fluctuates as the bead 
performs rotational Brownian motion. The FLWC theory accounts for these fluctuations via an 
effective boundary condition at the end(s) of the chain, which depends on applied force, bead 
radius, and the nature of the link joining the bead to the polymer (|l7l). We will study boundary 
conditions that are azimuthally symmetric; thus our end boundary conditions will be functions of 
t • z only. 

2.2 Fixed-angle bend 

We now suppose that our chain contains a permanent bend, whose location along the DNA, and 
angle, are fixed, independent of applied force. In this paper we will also neglect force-induced 
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unbinding of the deforming protein. (These effects are straightforward to incorporate into our 
analysis.) In addition, we neglect twist stiffness, which is legitimate since we wish to study a 
single bend in a polymer with unconstrained twist. (Twist stiffness effects will be important for 
experiments in which multiple bends occur or twist is constrained.) 

The Schrodinger-like equation Q must be modified by the inclusion of a "bend operator" 
which transforms tp at the bend. Suppose that the kink occurs at position s Q . Given the tangent- 
vector probability distribution ip at s Q — e (where e is infinitesimal), our goal is to determine 
^(t, s + e), the distribution just after the kink. If we denote the exterior angle of the kink by a, 
then t So _ (E -t So+e = cos a. Because twist is unconstrained, we may average over rotations; effectively, 
the bend occurs with uniform probability in the azimuthal angle: if t So _ e points directly along the 
z-axis, then t So+e is uniformly distributed in a cone at angle a to the z-axis. 

The bend operator then can be written using the kernel 

K a (i,i') = —5(i -t'-cosa). (6) 

ZTT 

The probability distribution ip(t', s — e) of tangent-vector angles just before the kink is related to 
the distribution ip(t, s Q + e) just after the kink by 

^(t, s + e) = J dl' K a (t, V) V(t', s - e). (7) 

Below (section l3.2j) we show that spherical harmonics diagonalize the operator ©. 



2.2.1 Distribution of bend angles 

Suppose that the bend occurs not for a single fixed angle, but a distribution of angles. We assume 
that t So _ e • t So+e = cos a = u is distributed according to the probability density function h(u), 
where h is normalized so that J dcf) J du h(u) = 1. Then the bend-operator kernel can be written 
as an integral over the probability distribution: 

If 1 

K h (t,t') = — du h(u) <5(t-t'-u). (8) 
2vr J_i 



3 Calculation 

The main quantity of interest in single-molecule experiments is the force-extension relation, which 
can be determined by solving equation Q for the tangent-vector probability distribution ip(t,s). 
The Schrodinger-like equation is solved using separation of variables in s and t, where the angular 
dependence is expanded in spherical harmonics (jlfT). 

oo 

^(t, S ) = ^%(s)^o(t). (9) 

3=0 

(By azimuthal symmetry, only the m = terms will enter in our formulae.) In the basis of spherical 
harmonics, the operator in equation Q is a symmetric tridiagonal matrix H with diagonal terms 



and off-diagonal terms 



n u = do) 



" + N /(2jTT)(2j+3) 
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The vector of coefficients at s is *&(s) = e sH ^(0) . UtI). This expression for ip(t,s) is exact if 
the infinite series of spherical harmonics is used. 

3.1 Force-extension relation 

Given the boundary conditions \&(s = 0) and \&(s = £), the partition function is 

Z = V T (s = e)e eH V(s = 0), (12) 
= Y,V j {s = £)[J B ) jk 9 h (a = 0). (13) 

The fractional extension of the chain is 

z IdlnZ 



L £ Of 



(14) 



We work in the ensemble relevant to most experiments, where the extension is determined for fixed 
applied force (different ensembles are not equivalent for single finite- length molecules (JzJ 2^, 2^)). 



Equation l)14j) applies for a chain of any length. However, we can show the structure of the partition 
function more clearly by sep arating InZ into two terms: one representing an infinite chain and a 
finite-length correction (|l7l). Let B = e H , denote by A* the largest eigenvalue of B, and let 
B = B/\*. Then B has eigenvalues with magnitude less than or equal to 1 and the logarithm of 
the partition function can be written 

In Z = £ln A* + ln[* T (s = £)B e V(s = 0)]. (15) 

Only the first term is considered in the usual WLC solution; the second term is the finite-length 
correction jl?! ) . Equation (|15|) is an exact expression for In Z which is difficult to evaluatae analyt- 
ically. We numerically calculate the force-extension relation by using equation (|15|) with the series 
truncated after N terms. This expression can be accurately numerically calculated, and the trun- 
cation error determined by comparing the results with different N. Our calculations use N = 30 
unless otherwise specified. 



3.1.1 Boundary conditions and bead rotational fluctuations 

The boundary conditions at s = and s = I affect the force-extension relation, because they 
alter the partition function as shown in equation (j!5j) . Some experiments appear to implement 
"half-constrained" boundary conditions, where the polymer is attached to a planar wall by a freely 
rotating attachment point, and the wall is perpendicular to the direction of applied force |2^). In 
this case the tangent vector at the end of the chain can point in any direction on the hemisphere 
outside the impenetrable surface (figure E[b)). The effects of different boundary conditions on the 
force-extension relation are considered in detail in reference (|l7h . In the "unconstrained" boundary 
condition the tangent vector at the end of the chain is free to point in any direction on the sphere 
(in 47T of solid angle, figure^). In this case ip(t) is independent of cos# and \& = (1, 0, • ■ ■ ,0). In 
the "half-constrained" boundary conditions (figure Hb), the tangent vector at the end of the chain 
can point in any direction on the hemisphere outside the impenetrable surface; then the leading 
coefficients of * are 1, 0.8660, 0, -0.3307, 0, 0.2073, 0. In the "normal" boundary condition, the 
tangent vector at the end of the chain is parallel to the z axis, normal to the surface (figure Eb). 
Then the coefficients of \l/ are all equal to 1 (j^). 

The FWLC formulation can also average over rotational fluctuations of spherical bead(s) at- 
tached to one or both ends of the polymer chain. The result is an effective boundary condition that 
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depends on applied force and bead radius (|17». Both the case of perpendicular wall attachment 
and bead attachment generate boundary conditions that are invariant under rotations about the z 
axis (the direction in which force is applied), and hence give boundary states of the form given in 
equation @. 

3.2 Bend operator 

We wish to represent the bend operator (equation (|7jl) in terms of spherical harmonics; the operator 
is diagonal in this basis. Denote x = t • t' and note that any function of x with —1 < x < 1 can be 
written as a series of Legendre polynomials |27l ): 

oo 

K a {x) = Y,hPi{x). (16) 

1=0 

The ki are determined by projecting the kernel K onto the Legendre polynomials, using the nor- 
malization relation J^ 1 Pi'(x)Pi(x)dx = ^rpr^z'- Therefore 

21 + 1 f 1 

ki = — j— J 5(x — cos a)Pi(x)dx, (17) 



47T 

21 + 1 



Pi(cosa). (18) 



4tt 

Next we use the addition theorem for spherical harmonics (|2 

A 1 

Pl{l ' V) = 2iTl ^ Y L(i')Yim(i), (19) 

m=—l 

Substituting equations (fTU|) and 1)18(1 in equation (fH)|) . we have 

oo I 

K a (i,i') = ^Pi(cosa) ^ YC m (i')Y lm (t). (20) 

Z=0 m=-Z 

Note that if a = 0, the kink operator reduces to the identity because Pi(l) = 1. 
The probability distribution tp just before the bend is 

oo j 

il>(t', *o - e) = E E *i^ s ° " (21) 
i=o fc=-i 

Note that in the case of azimuthal symmetry, the terms with k ^ are zero. To determine ip just 
after the bend, we substitute the expressions in equations (|20j) and (|21j) into the formula 

V>(t, So + e) = y dt' K Q (t, t') V(t', So - e). (22) 

The expression simplifies by the orthonormality of spherical harmonics: 

oo I OO j „ 

V(t, s + e) = P '( cos ") E E E "M s ° - £ ) / (23) 

i=0 m=-Z i=o fc=-j 

oo / OO j 

= ^Pi(cosa) ^ ilm(t)E E VjkiBo-eiSjiSmk ( 24 ) 

Z=0 m=-Z i=0 fc=-j 

oo Z 

= E E ^(cosa)#/ m (s -e)Yj m (t). (25) 

2=0 m=-Z 
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The transformation can thus be written ^fi m (s + e) = Pi(cos a)^i m (s — e). The probability 
distribution just after the kink differs from the distribution before the kink only in the multiplication 
of each term in the series by P;(cosa). We can represent the transformation by a diagonal matrix 
W such that 

*( So + e) = W*(s -e). (26) 

Because ip is azimuthally symmetric (only the m = terms appear in the series expansion), W has 
entries Wij = P/_i(cosa). 

3.2.1 Distribution of bend angles 

The representation of the bend operator in terms of spherical harmonics remains simple when 
the bend contains a distribution of angles described by h(u) (equation ©). As above, we expand 
Kh{x) in Legendre polynomials, Kh(x) = kiPi(x). The k\ are the projection of h(x) onto Legendre 
polynomials: 

h = -±- J dx h{x)P l {x). (27) 

The calculation is then identical to the case of a single bend angle, with the result ^i m (s + e) = 
ki^im(s — e). We can represent the transformation by a diagonal matrix Wh such that 

*( So + e) = W fc *( So -e). (28) 
3.3 Force-extension relation with bend 

Once the matrix W (which represents the bend operator in the basis of spherical harmonics) has 
been determined, calculation of the force-extension relation is straightforward. Suppose a single 
bend occurs at fractional position a along the chain. The partition function with a bend is 

Z b = * T (s = £)e {1 - a)m We aiH y{s = 0), (29) 

As before, we let B = e , denote by A* the largest eigenvalue of B, and define B = £?/A*. Using 
e iH = , the logarithm of the partition function is 

lnZ b = £ln\* + ln[* T (s = £) B^~ a)e W B a£ *(s = 0)]. (30) 

As above, the extension is z/L = l~ l din Z/df. 

4 Results 

Here we predict the magnitude of extension change induced by a single bend, in order to understand 
when such single-bending events will be experimentally detectable. We describe how the extension 
change induced by a bend depends on applied force, bend angle, contour length, and the position 
of the bend. 

In figure |31 we show the change in extension induced by a bend: the extension of the chain 
without the bend minus the extension of the chain with the bend. As expected, the extension 
change is larger when the bend angle is larger. In addition, we find that the change in extension 
has a maximum near an applied force of 0.1 pN. At this force, the change in extension due to the 
bend is a significant fraction of the persistence length (10-30 nm for A = 50 nm). 

[Figure 3 about here.] 
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For the largest bend angle, we show for comparison the prediction of Kulic et al. (Jl6|). The 
Kulic et al. result is valid in the high- force limit, and we find that their prediction and our result 
converge as the force becomes large. The Kulic et al. result is valuable because it is a simple 
analytical expression. Although our results are obtained numerically, they are valid over the entire 
force range. 

As the applied force increases, the polymer becomes more stretched and aligned with the force, 
decreasing the effect of the bend. For a classical elastic rod where thermal fluctuations are a weak 
perturbation the characteristic propagation length of elastic deformations is ^Jk^T A/F. Therefore, 
as the force increases, the region of the chain experiencing a significant deflection due to the bend 
drops. By this argument, one might expect that the largest change in extension due to the bend 
will occur for the lowest values of the applied force. However, as the force applied to the ends of the 
polymer goes to zero, the extension also approaches zero (on average, there will be no separation 
of the two ends). In this case the change in extension due to the bend approaches zero. The effect 
of the bend is therefore largest at intermediate force, where the molecule is extended by the force 
but not fully extended. 

We predicted the change in extension due to the bend with and without a bead attached to one 
end of the DNA, and for different values of the bead radius. In all cases, we predict similar values 
for the change in extension due to a bend (not shown). 

[Figure 4 about here.] 

In figure|l]we show how the change in extension induced by the bend varies with bend angle. The 
dependence of the extension change on angle is strong, suggesting that high-resolution experiments 
could measure the bend angle by measuring the change in extension due to a bend. For larger 
values of the applied force (F >1 pN), the result is independent of contour length of the polymer. 
However at low force (F=0.1 pN), where the change in extension due to a bend is largest, the 
results depend on the polymer contour length. 

[Figure 5 about here.] 

The dependence on the position of the bend is weak, unless the bend is within a few percent of 
one end of the polymer (figure EJ). We note that the curves in figure El are not reflection symmetric 
about the middle of the polymer. This occurs because we assume one end of the polymer (s = 0) is 
attached to a fixed surface, while the other end of the polymer (s = 1) is attached to a bead which 
can undergo rotational fluctuations. We chose to plot this case because it is a typical experimental 
geometry; in the case that both ends of the polymer experience identical boundary conditions, then 
the effects of a bend obey reflection symmetry about the middle of the polymer. 



5 Discussion 

We have described a theory of DNA elasticity applicable to bent DNA molecules. The finite worm- 
like chain model (FWLC) of polymer elasticity extends the WLC to polymers with L/A ~ 1 — 10 
(0). The FWLC includes chain-end boundary conditions and rotational fluctuations of a bead 
attached to the end of the polymer, modifications which are important for polymers with contour 
length a few times the persistence length. 

This work allows predictions of DNA force-extension behavior when a single bend occurs at a 
specified point along the chain. When the bend angle is constant (independent of applied force) the 
bend operator is diagonal in the basis of spherical harmonics, allowing straightforward calculation 
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of the effects of a bend. This mathematical description of a bend is suitable both for a bend with 
a single angle and for bends with a distribution of different bend angles. 

We demonstrate that the change in polymer end-to-end extension induced by the bend can be a 
significant fraction of the polymer persistence length: Az/A ~ 0.2 — 0.7 for bend angles of 90-180°, 
or Az ~ 10 — 35 nm for dsDNA, which has persistence length of approximately 50 nm. The change 
in extension due to the bend is predicted to show a maximum for applied force around 0.1 pN; 
for larger force the polymer conformation becomes highly extended and the influence of the bend 
decreases, while for low force the polymer extension approaches zero, independent of the presence 
of the bend. 

The alterations in polymer extension induced by the bend should be detectable in high-resolution 
single-molecule experiments. Since recent work in single-molecule optical trapping with DNA has 
demonstrated a resolution of a few nm(|2^. Eih. DNA extension changes of 10-35 nm due to a bend 
should be detectable. Furthermore, the predicted change in extension strongly depends on the 
bend angle, suggesting that high-resolution single-molecule experiments could directly estimate the 
angle of a protein-induced bend. 
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Figure 1: Typical experimental geometry of single- molecule force microscopy measurements. The 
DNA molecule is attached at one end to a surface and at the other end to a bead (radius R). A 
force F is applied to the bead, (a) DNA molecule in the absence of bound protein. The mean 
end-to-end extension is z. (b) DNA molecule with a single bend-inducing protein bound. The 
protein bends the DNA backbone through the external angle a at the bend site. As a result, the 
mean extension decreases to Zb en d- 
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Figure 3: Change in extension due to a bend as a function of force, determined by subtracting the 
extension of the chain with the bend from the extension without the bend. The contour length is 
(A) L=200 nm and (B) L=500 nm. For larger contour lengths, the prediction is similar to (B). 
We assume A = 50 nm, the bend is at the midpoint of the polymer, a bead of radius 250 nm is 
attached to one end of the chain, and half-constrained boundary conditions. For the largest bend 
angle, we show for comparison the prediction of Kulic et al. (16), which is valid in the high-force 
limit. 
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Figure 4: Change in extension due to a bend as a function of angle. (A) F=0.1 pN. (B) F=l pN. 
We assume A = 50 nm, the bend is at the midpoint of the polymer, a bead of radius 250 nm is 
attached to one end of the chain, and half-constrained boundary conditions. 
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Figure 5: Change in extension due to a bend as a function of the position of the bend along the 
polymer. (A) F=0.1 pN. (B) F=l pN. Note the difference in scale between the two panels. We 
assume A = 50 nm, L =500 nm, a bead of radius 250 nm is attached to one end of the chain, and 
half-constrained boundary conditions. 
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